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Abstract 

For a given density matrix p of a bipartite quantum system an 
asymptotical separability criterion is suggested. Using the continuous 
ensemble method, a sequence of separable density matrices is built 
which converges to p if and only if p is separable. The convergence 
speed is evaluated and for any given tolerance parameter k an iterative 
procedure is suggested which decides in finite number of steps if there 
exists a separable density matrix p K which differs from the matrix p 
by at most k. 

1 Introduction 

The notion of entanglement plays a central role in quantum communication 
and quantum computation. Usually the efforts are focused on quantifying 
entanglement itself, that is, describing the impossibility to prepare a state 
by means of LOCC (local operations and classical communications). One 
may, although, try to quantify separability rather than entanglement, this 
was shown to be a tool for geometrical classification of mixed multipartite 
entangled states 0. I dwell on the case of bipartite quantum systems. A 
state of such system is called separable if it can be prepared by LOCC. 
In terms of density matrices that means that p, its density matrix, can be 
represented as a mixture of pure product states. I suggest to replace finite 
sums of projectors by continuous distributions on the set of unit vectors. 

Recently continuous ensemble method was introduced |Hj which was 
applied to explore the separability of finite-dimensional quantum systems 
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HI EH- First, a geometrical characterization of robustly separable mixed 
states was provided in [§], then the separability problem was reduced to a 
finite number of numerical equations ^U]. This paper suggests an iterative 
procedure which, for any given tolerance parameter, decides if the state is 
separable or not for finite number of steps. 

This is done in the following way. In bipartite case the density operators 
are represented by distribution on the Cartesian product of unit spheres 
in subsystems' spaces. Given a density operator p, we consider it as an 
element of the space £ of all self adjoint operators in the product space 
25 = H®H. The density operator p is called robustly separable if it 
has a neighborhood of separable operators The robust separability of 
p is shown ^U] to be equivalent to the solvability of the following vector 
equation in £: 

VK = p 

with respect to the trace functional K on £ introduced in the first paper [H] 
of this series. When we fix a product basis in 23, this equation becomes a 
system of n transcendent equations with respect to n 4 variables 

In this paper an asymptotic procedure based on iterated function method 
is suggested which, in case if p is separable, approximates p by a sequence of 
separable density operators, otherwise this sequence diverges. The operator 
F is introduced 

F(X) = X + xfp-JJe^^lH/^l dSndS'^J 

where the integration is taken over the torus — the Cartesian product of unit 
spheres in H, H'. Then the sequence of operators Xq, X\, . . . , X m , ... is built 

' X = 

X x = F(X ) 
< ... 

X m = F(X m _i) 

i. * * * 

which converges if p is robustly separable. The constant A and the speed of 
this convergence are (roughly) evaluated, they depend on the dimension n 
of the state space and the tolerance k we choose. 
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2 Continuous optimal ensembles 



In this section I briefly summarize the basics of continuous ensembles ap- 
proach to quantifying quantum separability. In general context related to 
coherent states, continuous ensembles were considered yet by Glauber 
and Sudarshan Generalizing the fact that any convex combination of 
density operators is again a density operator, we represent density opera- 
tors as probability distributions on the unit sphere 1 in the state space 7i 
of the system (these ideas were also considered in 0). In particular, the 
well-known spectral decomposition of an operator is a special case of such 
representation, this will be addressed below. However the main feature of 
the suggested approach is that the distribution is taken over non- orthogonal 
states. 

Begin with the case of a single quantum system. Consider a probability 
distribution on C-B n whose density is a function /i(<j6), where <p £ < CB n ranges 
over all unit vectors in 7i. The density operator of this continuous ensemble 
is: 

p = J K<t>) \<f>M ds n (i) 

<t>&CB n 

where \4>)(4>\ is the projector onto the vector (<f>\ and dS n is the normalized 
measure on CB n . Effectively the operator integral p (^Q) can be calculated 
by its matrix elements. In any fixed basis {|ej}} in 7i, the matrix elements 
Pij = (ej| p \ej) are the following numerical integrals: 

Pij = (ei\p\ej) = J /i(0) (e* | <f>) (<j> | e,) dS n 

Smeared spectral decomposition. Since we are interested in robustly 
separable states, let us restrict ourselves to full-range density operators. The 
usual spectral decomposition p = J^Pk \ e k){ e k\ can be treated as an atomic 
measure on CB n whose density is the appropriate combinations of delta 
functions: 

x Pure states form a projective space rather than the unit sphere in TL. On the other 
hand, one may integrate over any probabilistic space. For technical reasons I prefer to 
represent ensembles of pure states by measures on unit vectors in TL. I use the Umegaki 
measure on CB n — the uniform measure with respect to the action of U(n) normalized so 
that J CB dS n — 1. Similarly, for bipartite case the integration will be carried out over 
the Cartesian product of unit spheres in appropriate state spaces. 
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^s Pec (0) = y^Pfc^o - e fc ) 

We introduce a 'smeared' version of the spectral decomposition. Let 
K be an integer number, then the density matrix p is represented as a 
continuous ensemble: 

P = ^Pfc|efc)(e fc | = J n{4>) \<j))(4>\ dS n 
The density \x was calculated in ^01, it equals 

Furthermore, the distribution (J2J) tends to the spectral decomposition of p 
as K tends to infinity. Denote by po the smallest eigenvalue of the density 
matrix p (recall that po > as we restrict ourselves to full-range density 
matrices) . Then the density \i in is positive for any such K that (K + 
l)n > 1/po- We may take 

K = — (3) 
npo 

Optimal entropy ensembles. It was shown in |5J HUj that, given a ro- 
bustly separable density matrix p in the product space OS = Ti ® Ti', there 
exists a self-adjoint operator X in 2? such that p can be represented as the 
following operator integral: 

P = J J e^ x ^ l^'X^'l dS n dS' n (4) 

</></>' ex 

Fix a product basis {(eje'j/]}- Then the condition @ takes the form 



( ei eV| P K-eV) = f f eW\xM) ( e . e >., | ^ ^ | g^V) dS n dS' n 

(5) 

and the separability problem reduces to finding the operator X from this 
system of equations. The essence of the separability problem is the existence 
of a solution of ©• In this paper I suggest an asymptotic procedure which 
gives the answer up to an arbitrary (but finite) precision. 
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3 Contraction mappings and the iterated function 
method 



A CONTRACTION MAPPING is any mapping F : V — > V in a metric space 
T> such that the distance d between any two points before the mapping is 
greater than the distance between them after the mapping: 

3C<1 Vx,yeV d(F(x),F(y)) < C- d(x,y) (6) 

The most important property of a contraction mapping is that there 
is exactly one point which is invariant under the mapping. The second 
crucial feature is that under sufficiently long iteration every point eventually 
contracts to this point. Using the definition ©, the speed of convergence 
can be evaluated from the value of the contraction parameter C and the 
diameter D of the space T>. After iV steps the distance between the fixed 
point of F and its iV-th approximation will not exceed C N D, therefore the 
given accuracy k will be achieved in at most 

In k — In D 



steps. A standard way to solve the equation f(x) 
function method is to define the function F(x): 

F(x) := x + X{a-f(x)) 

and choose a domain T> for x such that / acts within T> and a parameter 
A in order to make / contracting on T>. If the derivative /' is bounded 
on T>, then the contraction condition can always be satisfied by choosing 
sufficiently small A, namely, take any 



max | f 

In this case the contraction parameter in © is C = A/ max |/'|. 

Now let us apply this method to the equation © . Introduce the mapping 
F : £ -> £: 



(7) 

a using the iterated 



F(X) 



x + x 



\ 



(8) 



'el 
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with the domain T> = {X : < ^4} for a sufficiently large A, see section 
|I]for evaluations of the parameters in Q. If such X exists, then 

• we can always choose A which makes F contraction mapping (this is 
because \ x \^ > is bounded on V) 

• whatever be p, the zero operator is always in V 

This opens us the way to solve the separability problem. Start the it- 
eration procedure with the zero operator, then, if the solution exists, the 
sequence will converge, otherwise not. 

4 Evaluating the parameters 

In order to evaluate the parameters for the applicability of the iterated 
function method, we need to specify 

• The domain T> — to make the function F act within T> 

• The value of the parameter A — to make the mapping F contracting 

Assume that the density operator p is separable, then there exists a finite 
set of 1-dimensional product projectors {(f a f' a |} such that p is a convex 
combination 

P — ^ C-a |faf aXfaf a\ 

a 

The (real) dimension of the state of all self-adjoint operators in the prod- 
uct space 23 is n 4 . According to Karatheodori theorem, the projectors 
\faf'a){faf'a\ can be chosen so that their number in the above sum will 
not exceed n 8 . 

Now recall that we have a tolerance parameter k. What we need, is 
to assume that each product projector |f Q ,f' Q ,)(f Q ,f / a \ can be replaced by the 
mixture 
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Smeared spectral decomposition for the summands. Each sum- 
mand in @ is a product of density matrices in single-party spaces. Repre- 
sent it continuous mixture 



1^X01 dSn ® / |Ua(0') ^ (10) 



Let us first evaluate the range /^o < fJ- < Mi of the density for the single 
particle case. First note that the values of the uniform function f = ^2 c^r^ 1 
on the sphere ^ r\ = 1 ranges between 

</(ri,...,r n ) <c x 

where cq and ci are the least and the greatest coefficients, respectively. Then 
it follows from (J2J) that 



K-(2Tn)!n! n^™" 1 V (K + l)n / ~ K ■ {Kn)\ n\ \ (K + l)n y 
Denote 

((g + l)n)! 
Cx = K-(Kn)!n! 

For evaluation purposes replace pi — rR=qrn^ by 1. Substituting K = 1/npo 

from © and calculating p - ( K +i)n = n ' T^fe > n ' !i T a we obtain 

C K .^.^<»<C K (12) 

Recall that we consider our density operator up to a tolerance parameter k, 
hence we may assume from © that the smallest eigenvalue p$ in @ to be 
at least n/n 8 , therefore K can be evaluated from (j3J) as 

K = — = — (13) 

Continue the evaluation (j!2j) with po = k/u 8 : 

1 k 2 
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Then denote 

2n 14 n( ra8 /K) , . 

C A = - 2 (14) 

Taking into account that Ck, Ca >> 1, make the evaluation of /i rough but 
symmetric: 

< fi < C K ■ C A 



Ck ■ Ca 

then the magnitude of the density representing of each product p a in |(TU)) 
ranges within 



1 ^ 2 



CkCa 



< pt a (<f>) ■ A*a (<£') < C\C_ 



A 



Since the initial density matrix p is a convex combinations of the matrices 
p a , the matrix p itself can be represented as a continuous ensemble whose 
density ranges within the same bounds (|15|) . That means that the values of 
the operator X satisfy the following relations 

j_ y < e w\x\Ht) < C 2 kC 2 a (15) 

CkCaJ 

This gives us the domain for the function F we are going to iterate: 

||X|| < 2 In {C k Ca) (16) 

where the expression for Ck and Ca are obtained above in (fTT|) and 1)14(1 . 
respectively, and the norm is L\. 



Evaluating the derivative. In order to learn the value of the parameter 
A in (jHJ the derivative of F in any arbitrary direction must be evaluated. 
Let Y be a self-adjoint operator in 53, consider the operator F(X + tY) and 
calculate its derivative in the direction Y at t = 0: 



D 



d 
di 



F(X + tY) 



Y-X 



t=o 



,(<t>4,'\x\<fi<f>') 



Y 



'€% 



(17) 

Note that, according to 1(15)1. the values of e WI*W> are bounded by 
(CkCa) 2 , therefore, in order to evaluate the derivative, we take the maximal 
value, so the value of D is bounded 



dS n dS' 



S 



|£>ll < 



Y - A {C K C A f 



{^'\Y\4>cp') |#'X#'| ds„ds; 



'ex 



In order to evaluate the integral I use the formula obtained in |B] for the 
single particle case. For any self-adjoint operator A in TL 



A\ 



dSr 



<f>eCB n 



A + I-TrA 
n(n + 1) 



Recall that the L\ norm is used, this makes it possible to evaluate the traces 
of Y as they are bounded by the norm of Y, then 



A - J (<f>\A\ 

<f>£CB n 



dS n 



< \\A\\ ■ 1 



n(n + 1) 

Since the integral (|17|) is a linear function of the operator Y, 



(18) 



(#'|F|#') |#'X#'| dS n dS' n 



Y + I (8i Tri y + Tr n Y ® I + Tr Y ■ I ® I 
n 2 (n + l) 2 



•'eX 



where Tri,Trn stand for partial traces. So, taking into account (|TH)) 



l-DII < 



y - a (c K c A ) 



2 y + I ® y + Tr n y (8) I + Tr y • I g I 



n 2 (n + l) 2 



When we set A to be 



(CrCaY 



(19) 



we have the following evaluation of contraction constant \\D\\ < C\\Y\\ in 



C = 1 



n(n + 1) 



(20) 
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Evaluating the speed of convergence. According to (|16|) we have the 
following evaluation of the diameter of the domain X 



D < A\n[C A C K ) 

Then, fixing a tolerance parameter k and using the calculated contraction 
constant (|20[) we can apply ® to find the sufficient number of iterations 

= lnK-ln( 4ln(C A C K )) 
2m (l 



n(n+l) t 

and approximating the logarithm in the denominator we get 

N = 2n(n + 1) Hn(4 \n{C A C K )) + In ~ j (21) 

5 Concluding remarks 

Given a density matrix p in 23, a question arises if it is separable or not. 
When the dimension of at least one of spaces TL, TL' is 2, this question was 
given an effective and exact answer — the positive partial transpose (PPT) 
criterion due to Peres-Horodecki was suggested (HE*]. ^ ne criterion states 
that p is separable if and only if its partial transpose p Tn remains non- 
negative matrix. In higher dimensions PPT is only a necessary condition 
for a state to be factorizable as there exist entangled density matrices whose 
partial transpose if positive. 

In this paper the dimension of TL, TL' may take any finite value n. The 
separabilty problem is resolved with any given in advance finite precision. 
Let us summarize the procedure. 

• Given a density operator p in the product space 23 = TL (g) TL' ~ C n , 
calculate its minimal eigenvalue po. 

• Set up the tolerance parameter k and calculate the values of K (|13|) , 
C K (p} and C A d: 

n 7 ((K + l)n)\ 2n 14 n( n8 / R ) 

— ' ^ K ~ 1? — Tr? \ i — \ ' ^ A ~ 2 

k K ■ {Kn)\ n\ k z 

from which we derive the value of A (|19|) and the number of steps N 
HUD: 

A = 2 ; N = 2n(n + l) ( ln(Aln(C A C K )) + In -J 
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• Iterate the function F (JHJ) 



Xk+i 



F(X k ) 



( 

X k +\ p 

\ 




e 



starting from Xq = 0. At each step check ||AT||i and if it occurs that 
\\X\\i > 2 In (CkCa) that signals us that p is entangled. 

• After iV steps we get the approximation of p with accuracy n 

I emphasize that this method is not completely suited for practical cal- 
culations as the evaluations for the convergence parameters are very rough. 
This paper just demonstrates the principal possibility of using the continu- 
ous ensemble method to tackle the separability problem. 
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